Evidence for Asymptotic Safety from Lattice Quantum Gravity 
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We calculate the spectral dimension for nonperturbative quantum gravity denned via Euclidean 
dynamical triangulations. We find that it runs from a value of ~ 3/2 at short distance to ~ 4 at 
large distance scales, similar to results from causal dynamical triangulations. We argue that the 
short-distance value of 3/2 for the spectral dimension may resolve the tension between asymptotic 
safety and the holographic principle. 
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Introduction: The quantization of gravity is one of the 
great outstanding problems in theoretical physics. As 
is well known, a straightforward implementation of gen- 
eral relativity as a perturbative quantum field theory is 
not renormalizable by power counting, and the divergent 
counter-terms have been confirmed by explicit calcula- 
tion pp. Although it is still possible to formulate gravity 
as an effective field theory at low energies [2], at each 
order in the perturbative expansion new divergences ap- 
pear, requiring an infinite number of couplings that en- 
code the physics at high energy scales, leading to a loss 
of predictive power. 

Weinberg suggested in [3] that the effective quantum 
field theory of gravity might be asymptotically safe. In 
this scenario, the renormalization group flow of the cou- 
plings of the theory would have a nontrivial fixed point, 
with a finite dimensional ultraviolet critical surface of tra- 
jectories attracted to the fixed point at short distances. 
Thus, gravity would be ultraviolet complete and describ- 
able in terms of a finite, and possibly small, number of 
parameters, making it effectively renormalizable when 
formulated nonperturbatively. 

Evidence for the existence of a fixed point in the 4d the- 
ory of quantum gravity without matter has come mainly 
from lattice calculations [3H3 and continuum functional 
renormalization group methods [8lU3|. The renormal- 
ization group methods have shown the existence of a 
nontrivial fixed point in a variety of truncations of the 
renormalization group equations, indicating that the di- 
mension of the ultraviolet critical surface is finite. Var- 
ious truncations with more than three independent cou- 
plings suggest that the critical surface is three dimen- 
sional [T21 115] . 

Although the renormalization group studies are sug- 
gestive, the truncation of the effective action makes it 
difficult to systematically assess the reliability of the re- 
sults obtained using this method. A lattice formulation 
of gravity is thus desirable, given the possibility of per- 
forming calculations with controlled systematic errors. 
In a Euclidean lattice formulation the fixed point would 
appear as a second order critical point, the approach to 
which would define a continuum limit [14] . 



There exists an argument due to Banks [T5] (see also 
Shomer [TB]) against the possibility of the asymptotic 
safety scenario. The argument compares the density of 
states at high energies expected for a theory of gravity to 
that of a conformal field theory. Since a renormalizable 
quantum field theory is a perturbation of a conformal 
field theory by relevant operators, a renormalizable field 
theory must have the same high energy asymptotic den- 
sity of states as a conformal field theory. It follows from 
dimensional analysis, the extensive scaling of the quan- 
tities considered, and the fact that a finite-temperature 
conformal field theory has no dimensionful scales other 
than the temperature, that the entropy S and energy E 
scale as 

S~(RT) d -\ E~R d - l T d (1) 

where R is the radius of space-time under consideration, 
and T is the temperature. It follows that the entropy of 
a renormalizable theory must scale with energy as 

S~E^r. (2) 

For gravity, however, one expects that the high energy 
spectrum will be dominated by black holes. The d di- 
mensional Schwarzschild solution in asymptotically flat 
space-time has a black hole with event horizon of radius 
r d ~ 3 ~ GjyM, where Gn is Newton's constant and M 
is the mass of the black hole. The Bekenstein-Hawking 
entropy formula tells us that S ~ r d ~ 2 , so that 

S~E^. (3) 

This scaling disagrees with that of Eq. ([2| . Since consid- 
erable evidence supports the holographic argument lead- 
ing to Eq. p]), one is led to conclude that gravity can- 
not be formulated as a renormalizable quantum field the- 
ory. This is a potentially serious obstacle for asymptotic 
safety, and we revisit this issue below. 

Though they do not address this point, lattice calcula- 
tions with causal dynamical triangulations (CDT) have 
shown promising results [IHS] . In the parameter space of 
couplings a phase has been identified with a good semi- 
classical limit and a ground state consistent with (Eu- 
clidean) de Sitter space [B] [T7] . Also, there exist points 
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connected to this phase that are reasonable candidates 
for a second order critical point, though establishing the 
order of the transition at these points remains a difficult 
numerical problem [THl HH] . 

The original formulation of Euclidean dynamical trian- 
gulations (EDT) is a variant of quantum Regge calculus 
[2"0] where the lattice geometries entering the Euclidean 
path integral are approximated by triangles of fixed edge 
lengths. The dynamics is contained in the connectiv- 
ity of the triangles and their higher dimensional analogs, 
called simplices. The dimension of the simplices is fixed, 
but taking four-dimensional building blocks does not en- 
sure a four-dimensional geometry, and the effective frac- 
tal dimension must be determined from the simulations. 
There are two phases in the simplest implementation of 
the model. Unfortunately, neither of these phases resem- 
bles semiclassical general relativity in four dimensions. 
Also, the critical point separating these two phases was 
shown to be first-order, thus ruling out the possibility of 
taking a continuum limit [211 E2] . 

The CDT approach was introduced by Ambjorn and 
Loll mainly in response to these difficulties [23 . The 
good results obtained so far with the CDT method have 
been attributed to the introduction of a causality condi- 
tion, where one distinguishes between spacelike and time- 
like links on the lattice so that an explicit foliation of the 
lattice into spacelike hypersurfaces can be introduced, 
where all the hypersurfaces have the same topology. Only 
geometries admitting such a global foliation are included 
in the space of triangulations defining the measure of the 
path integral. In addition to the emergence of de Sit- 
ter space [5J, the CDT calculations lead to the striking 
result that the effective fractal dimension (the spectral 
dimension) runs as a function of distance [24 , in accord 
with renormalization group calculations that support the 
asymptotic safety scenario [25] . 

Inspired by these successes, we revisit the original EDT 
formulation in this work. Our motivations are the follow- 
ing: First, the fixing of the foliation in CDT is potentially 
at odds with general covariance. If consistent with gen- 
eral covariance, it should be possible to obtain the same 
results using an EDT formulation, which is explicitly co- 
variant from the outset. 

Second, there are three parameters that enter the bare 
lattice action in CDT, whereas the original EDT studies 
typically included only two. The renormalization group 
results suggest that the ultraviolet critical surface is three 
dimensional, so we should consider the possibility that 
the key ingredient of the CDT formulation is the third 
coupling, rather than the causality condition. Both EDT 
and CDT have the bare Newton and cosmological con- 
stants. CDT also introduces an anisotropy ratio of the 
lengths of timelike and spacelike links, so that CDT ge- 
ometries are anisotropic along the time direction [5J . Our 
aim is to study how EDT behaves when a third coupling 
is introduced. 



Method: Within the EDT formalism, the path inte- 
gral for Euclidean gravity becomes the discrete partition 
function [20J [26] 




where Ct is a symmetry factor that divides out the num- 
ber of equivalent ways of labeling the vertices in the tri- 
angulation T. The Euclidean Einstein-Regge action is 

S E = -k 2 N 2 + k 4 A 4 , (5) 

where iVj is the number of simplices of dimension i, and 
where n 2 and K4 can be related to the bare Newton's 
constant Gn and the bare cosmological constant A. The 
term in brackets in Eq. Q is a nontrivial measure term 
[27] . where the product is over all two-simplices (trian- 
gles), and o(tj) is the order of triangle j, i.e. the number 
of four-simplices to which the triangle belongs. We vary 
the free parameter j3 as an additional independent cou- 
pling constant in the action, bringing the total number 
of couplings to three. Most of the previous work on EDT 
considered the partition function in Eq. Q with f) = 
only. 

Almost all of the work on four-dimensional EDT used 
geometries constructed of combinatorial triangulations, 
in which each d-simplex is defined by a unique set of 
d+1 distinct vertices, and each d-simplex has a set of 
unique neighbors. It is possible to relax this constraint, 
so that distinct simplices can be defined by the same set 
of vertices and the neighbors of a d-simplex are no longer 
unique [28] . It was also shown in Ref . [28] that one must 
keep the restriction that each four-simplex is defined by 
a set of unique vertices. These restrictions define a set 
of degenerate triangulations. It has been shown numeri- 
cally that using degenerate triangulations leads to a sig- 
nificant reduction in finite-size effects (~ a factor of 10) 
for both two-dimensional EDT models, where analytical 
results are available for comparison |14j . and for four- 
dimensional models in the unphysical phases where EDT 
has been studied extensively using combinatorial triangu- 
lations [5S]. We use the set of degenerate triangulations 
in the sum of Eq. Q, though we also use simulations 
with combinatorial triangulations as a cross-check. 

The numerical implementation of EDT is well- 
documented, and we refer elsewhere for more details [14] . 
In brief, we implement a Monte Carlo integration over 
a set of four-dimensional degenerate triangulations with 
the partition function of Eq. Q and fixed topology S A . 
We use the Metropolis algorithm with a set of local up- 
date moves. It is convenient to keep the four- volume 
approximately fixed, so we include a term in the action 
8 A I N{ — A4I to keep the four- volume close to the value 
n[. The bare cosmological constant (or equivalently K4) 
is tuned to its critical value, so that one can take an 
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FIG. 1: Schematic of the phase diagram as a function of K2 
and /?. 

infinite- volume limit. This leaves a two dimensional pa- 
rameter space, which is explored by varying k,i and ft. 
Our code is new, but we have checked that we are able 
to reproduce many results from the extensive literature 
on combinatorial triangulations and the work of Ref. [25] 
for degenerate triangulations. 

A schematic illustration of the phase diagram for EDT 
is shown in Fig. [T] within the k 2 , ft plane. The line at 
ft = has received much attention in previous studies, 
with a first-order phase transition separating the col- 
lapsed and branched polymer phases being fairly well- 
established. Previous work 26J suggests that adding a 
nontrivial measure term with sufficiently negative j3 leads 
to the appearance of a third phase, which we label "ex- 
tended" in Fig. [T] Our simulations with degenerate and 
combinatorial triangulations thus far appear to support 
this picture. Although it is nearly certain that line AB is 
first order, it is not yet clear from our simulations what 
the order of the transition is along the dotted lines BC 
and BD in Fig. [T] The dotted lines indicate that these 
transitions appear to be softer than the transition across 
AB, and it is possible that one or both of these may be 
crossovers. If this is the case, then the end point of line 
AB would still be a candidate for a second order critical 
point. Determining the order of such phase transitions is 
difficult, requiring rather large lattices, and the diagram 
in Fig. [TJ should be regarded as tentative. As a first step, 
we have begun to characterize the extended region to see 
if it has the behavior expected of a four-dimensional the- 
ory of gravity, postponing a detailed study of the critical 
behavior and the existence of a continuum limit. 

To determine whether the extended phase has semi- 
classical gravitational properties, we compute the spec- 
tral dimension. We follow the approach used for CDT in 
Ref. [23] . The geometries that dominate the path-integral 
are presumably not the smooth manifolds of classical 
gravity, but rather they are like the paths that dominate 
the path-integral in nonrelativistic quantum mechanics: 
continuous, piecewise linear geometries that are nowhere- 
differentiable and possess a fractal structure. The spec- 
tral dimension defines the effective dimension of such a 
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FIG. 2: Spectral dimension as a function of diffusion time a, 
including a fit to Eq. (JtJ . The width of the band shows the 
statistical error in the fit. 

fractal geometry via a diffusion process, and it reduces to 
the usual definition of dimension on smooth manifolds. 

The spectral dimension of the lattice geometries of the 
type considered in this work can be obtained from |24| 



, <flog(P(<r)) 
(floger 



(6) 



where (P(a)) is the ensemble average of the return proba- 
bility, i.e. the probability that a random walk will return 
to its origin after a steps. This formula receives finite- 
size corrections when a becomes much larger than the 
lattice volume N^ Ds . The spectral dimension is not the 
only way to define an effective dimension on a fractal 
and in general will differ from other definitions of di- 
mension, such as the Hausdorff dimension. The spectral 
dimension warrants particular attention because it is the 
relevant scaling dimension entering the thermodynamic 
equations of state on a fractal geometry [2H] . 

Results: Figure [2] shows our new result for the spectral 
dimension plotted as a function of a. This is the cen- 
tral result of our letter. We obtained this result from an 
ensemble of degenerate triangulations with a volume of 
4000 four-simplices in the extended phase with ni = 2.1 
and = —1.0. We analyzed ~ 1000 configurations us- 
ing the average of the return probability as in Eq. (§), 
and the errors are determined from a single elimination 
jackknife procedure. To reduce statistical errors the data 
were binned in bin sizes of two from a = 2 to 80 and 
in bin sizes of 4 from u = 80 to 288. The fit is to the 
functional form suggested by Ref. |24) . 



(7) 
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with the constants a, b, and c determined by the fit. 

Our preferred fit is shown in Fig. [2] with statistical 
errors. This fit uses the full covariance matrix in the 
estimate of \ 2 with a ranging from 10 to 146 in steps 
of 4, yielding a x 2 /dof = 35/32 and a confidence level 
(corrected for finite sample size) of CL= 0.37. Variations 
of the fit range and fit function were used to estimate a 
systematic error on the asymptotic value of Ds, always 
under the assumption that Dg approaches a constant as 
a goes to infinity. The results of our preferred fit are 
D s (oo) = 4.04 ± 0.26, D S (Q) = 1.457 ± 0.064, where the 
errors include both the statistical error and a systematic 
error associated with varying the fit function and the fit 
range added in quadrature. As one can see from the 
small a data, -Ds(O) is not very sensitive to the choice 
of fit function, nor to changing the minimum value of a 
to values up to 40. The behavior of the spectral dimen- 
sion is generic within the extended phase and does not 
depend on the precise values of the bare parameters used 
in the action. Finally, we remark that calculations of 
the spectral dimension with combinatorial triangulations 
in the extended phase give similar results, but require 
significantly larger lattice volumes. 

Conclusions: Our results for the running spectral di- 
mension are in reasonable agreement with the CDT re- 
sults: D s {oo) = 4.02 ± 0.1 and D s (0) = 1.80 ± 0.25 gj. 
This supports the idea that the explicitly covariant EDT 
formulation is a valid candidate for a theory of quantum 
gravity, with a spectral dimension ~ 4 at long distances, 
and it provides nontrivial evidence that CDT and EDT 
may be in the same universality class. 

Our result for the short-distance spectral dimension 
Ds(0) = 1.457 ± 0.064 is inconsistent with the renor- 
malization group result that Ds{0) — 2 exactly [2"S"j . 
However, returning to the holographic argument against 
asymptotic safety, we note that for the asymptotic den- 
sity of states one should use the short-distance spectral 
dimension as the scaling dimension in Eq. ^. Under 
the plausible assumption that the relevant dimension in 
the holographic scaling argument is also the spectral di- 
mension, we note that Eqs. |2]) and Q agree precisely 
when Ds = 3/2. Then, our result from EDT resolves 
the tension between asymptotic safety and holographic 
entropy scaling. In light of this, it is worth revisiting the 
renormalization group arguments leading to -Ds(O) = 2, 
as well as embarking on more precise CDT calculations 
to see whether they confirm the value 3/2. If one could 
demonstrate the existence of a continuum limit at a sec- 
ond order critical point, this result could be the first clue 
as to how a renormalizable quantum field theory of grav- 
ity may yet be consistent with the holographic principle. 
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